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ABSTRACT 


Atmospheric  turbulence  introduces  random  phase  distortions  in  optical 
imaging  systems.  The  development  of  new  laser  and  imaging  systems  requires 
information  on  the  spatial  and  temporal  distribution  of  this  atmospheric 
turbulence.  Measurements  of  the  image  spread  and  the  jitter  induced  by  the 
atmosphere  on  an  optical  system  provide  two  techniques  to  quantify  these 
phenomena.  This  thesis  evaluates  a  Spectra  Sources  Lynxx  PC  Plus  charge 
coupled  device  (CCD)  array  as  an  atmospheric  turbulence  sensor.  Data 
acquisition  and  processing  programs  were  written  to  measure  the  image  spread 
of  a  point  source  and  centroid  jitter  of  a  point  source  imaged  through  the 
atmosphere.  Since  atmospheric  jitter  measurements  require  high  image  frame 
rates  on  the  order  of  200  images  per  second,  a  large  portion  of  this  thesis 
involved  measurements  of  the  times  for  the  CCD  detector,  interface  board  and 
IBM  compatible  computer  to  perform  their  tasks.  Recommendations  for  higher 
performance  are  presented. 
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I  INTRODUCTION 


This  thesis  considers  a  method  used  to  measure  atmospheric  turbulence 
using  a  silicon  charge  coupled  device  (CCD)  camera.  Light  from  a  distant  star 
consists  of  almost  perfect  plane  waves  before  entering  the  earth’s  atmosphere. 
When  the  waves  traverse  the  atmosphere,  local  variations  in  the  index  of 
refraction  induce  warping  of  the  wave  fronts.  Therefore,  a  telescope  pointed  at 
a  star  will  produce  a  distorted  image,  whose  spatial  and  temporal  characteristics 
depend  upon  the  properties  of  the  atmospheric  path  followed  by  the  starlight. 
Typically,  the  atmosphere  introduces  both  image  spread  and  temporal  jitter,  so 
that  a  CCD  array  placed  at  the  telescope’s  focal  plane  will  reveal  a  highly 
unstable  and  complex  irradiance  pattern.  It  has  been  found  that  both  a  high 
frame  rate  (>  200  per  second)  and  correspondingly  short  exposure  time  (<  5 
ms)  are  necessary  to  avoid  under-sampling  the  dynamic  variations  of  typical 
images. 

The  goal  of  this  thesis  was  to  find  a  way  to  measure  atmospheric  turbulence 
with  a  CCD  camera  by  optimizing  the  speed  and  efficiency  of  the  camera 
system’s  software  and  hardware.  The  software  was  written  in  the  C  language. 
Software  components  that  can  affect  system  speed  include  both  the  various 
algorithms  used  and  the  overall  efficiency  of  the  compiled  program  code.  Three 
compilers’  outputs  were  compared  for  executable  code  speed  and  utility: 
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Microsoft  Quick  C  2.5,  Borland  C++  2.0  and  Borland  Turbo  C++  1.0.  Fixed 
hardware  components  included  a  Texas  Instruments  TC21 1  CCD  image  sensor 
in  a  Spectra  Source  PC  Lynxx  Plus  camera  attached  to  an  8  bit  IBM  PC 
compatible  interface  card.  The  camera  system  was  installed  in  two  IBM 
compatible  personal  computers,  one  using  an  Intel  80386  and  the  other  a  80486 
processor,  and  the  overall  speeds  were  compared.  Finally,  the  camera-computer 
interface  and  the  CCD  image  sensor  were  analyzed  to  determine  whether  such 
a  low-cost  commercial  system  would  be  suitable  for  turbulence  measurements, 
or  whether  a  customized  device  would  have  to  be  fabricated. 


* 
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II  BACKGROUND 


A.  ATMOSPHERIC  TURBULENCE 

Turbulence  in  the  atmosphere  causes  many  problems.  In  astronomy,  it 
introduces  distortion  that  obscures  image  detail.  Since  exoatmospheric  use  is 
not  practical  for  most  telescopes,  knowing  the  spatial  and  temporal  distribution 
of  turbulence  is  essential  in  the  planning  stages  of  a  new  facility  or  in  diagnostic 
evaluation  of  test  systems.  Figure  1  sketches  a  wave  front  incident  on  a  turbulent 
region  [Ref.  1],  Random  phase  fluctuations  in  the  index  of  refraction  field 
produce  scintillations  and  image  blurring.  Twinkling  or  scintillation  arises  from 
the  interference  of  starlight  that  traverses  multiple  paths  through  the  atmosphere. 
Image  blurring  arises  from  the  reduction  in  spatial  coherence  of  the  wave  front. 


B.  PARAMETERS  FOR  THE  MEASUREMENT  OF  TURBULENCE 

The  atmospheric  index  of  refraction  depends  on  both  temperature  and 
pressure.  Since  pressure  fluctuations  disperse  rapidly,  at  the  speed  of  sound, 
temperature  fluctuations  are  the  main  contributor  to  atmospheric  optical 
turbulence.  Structure  functions  provide  a  way  to  quantify  the  statistics  of 
atmospheric  turbulence.  For  propagation  through  the  atmosphere,  optical 
parameters  will  involve  an  integral  of  the  refractive  structure  function  along  the 
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Figure  1 .  Optical  atmospheric  turbulence  distorts  a  coherent  phase  front, 
optical  path.  The  primary  optical  parameter  is  the  spatial  coherence  length,  r0. 

It  represents  the  transverse  autocorrelation  length  of  the  electromagnetic  field. 

A  related  parameter  is  the  isoplanatic  angle,  0o.  It  represents  the  angle  between 
two  different  paths  where  the  Strehl  ratio  of  an  ideal  adaptive  optics  system  is 
within  e1  of  perfect  correction  (Ref.  2].  Another  parameter,  the  Greenwood 
frequency,  fg,  represents  the  electrical  bandwidth  needed  to  remove  atmospheric 
phase  distortions  with  an  adaptive  optical  system.  * 
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1.  Structure  Function 


Atmospheric  turbulence  produces  localized  variations  in  the  index  of 
refraction  along  an  optical  path.  Tatarski  [Ref.  3]  defines  the  structure  function 
by 


D*<r2~ri)  =  Ux(r2)  ]  2>,  (l) 

where  <  >  represents  an  ensemble  spatial  average,  7,  and  7Z  represent  the 
location  of  two  points  in  space,  and  x  represents  an  atmospheric  parameter  of 
interest,  such  as  temperature  or  index  of  refraction.  Assuming  the  turbulence  is 
isotropic,  homogeneous  and  incompressible  [Ref.  4]  the  Kolmogorov  theory  of 
turbulence  [Ref.  5]  shows  that 


Dx  =  cUU\  (2) 

where  Cx2  is  a  structure  parameter  and  r12  is  the  distance  between  the  two  points 
in  space.  By  combining  equations  (1)  and  (2),  the  thermal  structure  parameter, 
CT2,  which  characterizes  the  mean  squared  temperature  difference  between  two 
points  in  space  is  [Ref.  3], 


c;  = 


((Tj-Tj)2) 


-2/3 

r12 


(3) 


where  T2-T,  is  the  temperature  difference  and  r12  is  the  distance  between  the  two 
points  in  space.  There  is  a  similar  expression  for  the  index  of  refraction  structure 
parameter,  Cn2.  Taking  the  partial  derivatives  of  the  atmospheric  index  of 
refraction  [Ref.  3], 
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n-1 


79x10  ~6P 
T 


(4) 


where  P  is  the  atmospheric  pressure  in  mbar  and  T  is  the  Kelvin  temperature, 
and  assuming  isobaric  turbulence  allows  us  to  write  Cn2  in  terms  of  Cx2  [Ref.  6], 

d  =  (  Y"  —  )  2  c\.  (5) 

*T»2 

The  pressure  is  assumed  to  be  constant,  since  small  random  pressure 
differences  disperse  rapidly. 

2.  Spatial  Coherence  Length 

The  spatial  coherence  length,  r0,  as  stated  above,  measures  the 
transverse  autocorrelation  length  of  a  wave  front.  It  represents  the  effective 
aperture  diameter  of  a  diffraction-limited  optical  system  with  a  similar  angular 
resolution  as  the  system  under  study.  Typical  values  of  rD  range  from  two  or 
three  centimeters  for  high  turbulence  to  as  much  as  thirty  centimeters  for  low 
turbulence,  when  measured  from  the  earth’s  surface  upward.  When  observing 
a  star  from  the  ground  with  an  optical  device,  the  Cn2  component  of  r0  is 
cumulative  along  the  optical  path,  even  though  most  of  the  degradation  of  the 
wave  front  occurs  close  to  the  ground.  Integrating  the  optical  turbulence  Cn2(z) 
along  the  path  [Ref.  6], 

ra  =  [ k2  sec$  cj;  U)  Wz)  dz]  '3/5,  (6) 

where  k  is  the  wave  number(2rr/X),  $  is  the  zenith  angle,  L  is  the  vertical  path 
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length  and  W(z)  is  a  weighting  function.  Typically  W(z)  has  one  of  three  forms: 
for  a  plane  wave,  W=1;  for  a  spherical  diverging  wave,  W=(z/L)s/3;  and  for  a 
converging  wave,  W=  (1 -z/L)5/3.  The  factor  of  sec  0  compensates  for  the  slant 
path  through  a  horizontally  stratified  atmosphere. 

Temperature  fluctuations  along  a  vertical  path  have  been  measured  by 
two  means:  by  microthermal  probes  carried  by  a  meteorological  balloon  and  by 
an  echo  sounder  [Ref.  7  and  8].  Both  of  these  methods  provide  credible  values 
of  r0.  However,  there  is  a  more  direct  method  for  calculating  the  coherence 
length  at  the  earth’s  surface  using  an  optical  point  source  irradiance  distribution 
[Ref.  9],  The  discrete  irradiance  values  from  a  two  dimensional  irradiance 
distribution  represent  the  point  spread  function,  P(x,y).  The  line  spread  function 
for  x,  L(x),  is  the  summation  of  all  of  the  y  irradiance  values  from  P(x,y)  for  each 
x.  L(y)  is  a  similar  summation  of  the  x  values  for  each  y.  Then  take  the  Fourier 
transform  of  L(x)  or  L(y)  to  get  the  optical  transfer  function  (OTF).  Assuming  the 
turbulence  is  isotropic,  that  is  symmetric  with  respect  to  rotation  about  the  path, 
and  that  it  is  laterally  stationary,  the  absolute  value  of  the  OTF  equals  the 
modulation  transfer  function  (MTF).  For  an  optical  system  imaging  a  star  the 
observed  quantity,  MTF0,  is  a  product  of  the  MTF’s  of  the  source  or  star,  the 
atmosphere  and  the  instrument  [Ref.  9], 

MTF0  =  MTFS  'MTFa  'MTFi  .  (7) 

MTFt  for  a  star  corresponds  to  a  point  source  at  infinity  and  is  unity.  MTFj  is  a 
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measured  instrumental  transfer  function.  To  determine  the  atmospheric 
component,  MTF,,  from  the  observed  quantity,  MTFe,  we  divide  MTF0  by  the 
measured  instrumental  transfer  function,  MTFj. 

Fried  shows  that  atmospheric  MTF,  is  a  function  of  the  atmospheric 
wave  structure  function,  which  is  an  integral  of  the  index  of  refraction  structure 
parameter  along  the  path  [Ref  2].  MTF,  reduces  to  a  form  [Ref.  9], 

-3.44  ( 

MTFa  =  e  r°  (8) 

where  R  is  the  effective  focal  length,  X  is  the  wavelength  of  light  and  v  is  the 
spatial  frequency.  To  determine  the  atmospheric  parameter  rOI  we  took  a  one 
dimensional  fast  Fourier  transform  of  a  line  spread  function  of  a  star  image. 
Dividing  this  by  the  known  instrument  function  gives  the  spatial  spectrum  similar 
to  Equation  (8).  The  next  step  was  to  determine  the  e'1  frequency  for  this 
experimental  spatial  spectrum  from  which  [Ref.  9], 

rQ  =  2  .  l  X  R  v  e,  (9) 

where  ve  was  the  e 1  frequency  of  the  MTF,. 

Yet  another  way  to  calculate  r0  uses  image  centroid  jitter,  also  called 
beam  wander,  by  calculating  the  x  and  y  centroid  standard  deviation  from  a 
sequence  of  point  source  images.  It  is  equal  to  aF,  where  a  is  the  atmospheric 
angle  of  arrival  fluctuations  and  F  is  the  optical  system’s  effective  focal  length. 
The  angle  of  arrival  is  sketched  in  Figure  2.  Mathematically,  angle  of  arrival 
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Pt*n* 

Wave 

Figure  2.  The  angle  of  arrival  is  the  least-squares  plane  that  fits  the  electric  field 
vector. 


fluctuations  are  equivalent  to  least-squares  planes  that  fit  the  electric  field  across 
an  aperture.  In  terms  of  the  index  of  refraction  structure  parameter  [Ref.  5  and 
10], 


«2  =  1-°>23-  /LC2(z)  dz,  (10) 

D1/3  Jo 

where  D  is  the  diameter  of  the  aperture  and  L  is  the  optical  path.  Combining 
equation  (10)  with  equation  (6)  for  a  plane  wave, 

re  =  [1 . 418lc2sec<t>  D1/3a2]  '3/5.  (11) 

To  calculate  r0  from  the  centroid  motion  requires  a  very  high  image  frame  rate  to 
"freeze"  the  atmosphere.  Characteristically  at  least  200  frames  per  second  are 
required  in  order  to  avoid  aliasing  that  would  otherwise  underestimate  the  true 
centroid  variance. 
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3.  Isoplanatic  Angle 

The  isoplanatic  angle  measures  the  angular  coherence  in  the  vicinity 
of  an  object.  It  defines  a  cone  that  constrains  the  angles  over  which  an  adaptive 
optical  system  will  provide  valid  correction.  In  terms  of  Cn2  [Ref.  1], 

e0  =  [2.91  k2  flcZ(z)  z5l3dz]-3/s,  (12) 

Jo 

where  the  parameters  are  the  same  as  those  of  r0.  Simple  instruments  exist  to 
measure  0O  [Ref.  11].  For  this  reason,  only  the  coherence  length  is  addressed 
in  this  thesis. 

4.  Greenwood  Frequency 

The  Greenwood  frequency  was  introduced  in  an  earlier  paper  by  Darryl 
Greenwood  [Ref.  12  and  13],  It  is  a  mean  temporal  frequency  of  an  adaptive 
optics  system  and  depends  on  the  wind  velocity,  V,  and  Cn2  along  the  path.  For 
a  Kolmogorov  turbulence  spectrum,  the  Greenwood  frequency  is  [Ref.  13], 

f  =  2.31A'6/5  [JLC$(z)  V5/*(z)  dz] 3/5 ,  (13) 

where  k  is  the  wavelength  of  light,  z  is  the  distance  along  the  propagation  path, 
V(z)  is  the  wind  speed  and  L  is  the  distance  from  the  source  to  the  receiving 
optical  device. 
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C.  CHARGE  COUPLED  DEVICE  (CCD) 


An  atmospheric  coherence  length  sensor  needs  an  imaging  detector  to 

t  :j 

measure  the  point  spread  function.  The  speed,  reliability  and  availability  of 
charge  coupled  devices  (CCD's)  makes  them  worthy  of  consideration.  The  basis 
of  a  CCD  is  a  metal-oxide-semiconductor  (MOS)  capacitor  forming  potential  wells 
and  channels  that  comprise  light-sensitive  pixels  and  read-out  registers.  Figure 
3  is  a  sketch  of  the  MOS  capacitor  structure  [Ref.  14].  Electrons  excited  by  the 
photoconductive  effect  are  trapped  in  potential  wells  formed  under  the  positively 


Figure  3.  The  metal-oxide-semiconductor  (MOS)  capacitor  structure, the  basis  of 
the  charge  coupled  device  (CCD). 

charged  metal  gate  contact.  The  number  of  electrons  trapped  is  proportional  to 
the  integrated  irradiance  during  an  exposure. 
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There  are  several  methods  to  read-out  the  device.  Voltages  with  different 
phases  placed  on  the  gates  transfer  charge  from  pixel  to  pixel  or  potential  well 

-j 

to  potential  well.  A  set  of  electrodes  connected  together  is  called  a  phase. 
Figure  4  shows  an  example  of  a  three  phase  device  with  three  pairs  of  gates  (G) 
and  three  lines  (L)  connecting  the  electrodes  [Ref.  14].  A  positive  electrical  pulse 


Figure  4.  Three  pairs  of  gates  (G)  are  connected  with  the 
lines  (L)  to  form  three  different  phases. 

from  a  clock  forces  a  transfer  of  charge  from  one  phase  of  pixels  to  the  next. 
Clocking  a  phase  of  cells  reduces  the  barrier  between  pixel  columns  allowing 
charge  transfer.  A  bucket  of  water  running  down  several  steps  with  boards  or 
barriers  on  each  step  is  a  visual  representation  of  this  process.  The  water  flows 
down  to  a  lower  potential  each  time  the  board  is  lifted  up.  The  number  of 
phases  used  to  transfer  out  charge  packets  represents  the  number  of  steps. 

There  are  two  methods  used  to  insure  that  efficient  charge  transfer  occurs 
in  one  direction.  The  first  uses  multiple  phases  to  separate  the  charge  packets. 
The  Texas  Instruments  TC21 1  [Ref.  15]  uses  the  second  method  which  is  to  force 
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an  asymmetry  in  each  well  using  an  ion  implantation  zone  between  each  pixel. 
The  clocking  method  is  the  monophase  mode,  also  referred  to  as  a  1  +  Vs  phase 
mode,  sketched  in  Figure  5  [Ref.  16].  One  phase  has  an  intermediate  voltage 

r 

level  while  the  voltages  applied  to  the  other  phase  vary  on  both  sides  of  this 
level. 

The  charge  packets  transfer  down  channels.  There  are  two  types  of 
channels,  surface  and  buried.  The  water  analogy  also  works  here.  In  a  buried 
channel,  there  is  less  loss  of  charges,  similar  to  a  pipe.  A  surface  channel, 
however,  has  a  higher  capacity,  like  a  canal,  but  this  CCD  has  higher  noise  from 
surface  imperfections. 


Figure  5.  The  monophase  clocking  method  holds  a  phase  at  an  intermediate 
level  while  the  voltages  are  applied  on  either  side. 
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During  the  read-out  of  a  CCD,  image  smear  will  occur  since  light  will  still 
produce  charge  in  the  moving  CCD  rows,  unless  a  shutter  is  incorporated. 
Different  techniques  exist  to  reduce  image  smear  during  read-out.  Frame  transfer 
devices  overcome  this  disadvantage  by  effectively  having  two  arrays,  one  for 
image  exposure  and  a  second  for  storage.  Moving  charges  from  the  image  area 
to  an  opaque,  masked  set  of  pixels  quickly  allows  another  exposure  to 
commence  while  reading  the  storage  area.  Figure  6  shows  two  examples  [Ref. 
14]. 


J 
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Figure  6.  Two  schemes  for  CCD  read-out.  Top,  the  interline  method;  bottom, 
frame  transfer  method. 


15 


Ill  DISCUSSION 


A.  EQUIPMENT 

This  thesis  investigated  a  Spectra  Source  Lynxx  PC  Plus  CCD  Imaging 
System  [Ref.  17  and  18].  The  original  intent  was  to  use  this  CCD  camera  with 
a  telescope  and  personal  computer  to  develop  software  to  measure  the 
atmospheric  coherence  length,  r0.  During  this  process  it  became  clear  that  the 
off-the-shelf  CCD  system  could  not  provide  sufficient  frame  rates  to  achieve  the 
intended  goal.  Consequently,  the  task  shifted  to  determine  the  causes  of  the 
sluggish  acquisition  rates.  The  CCD  detector  array,  computer  interface  card  and 
the  IBM  compatible  PC  were  investigated. 

1.  CCD  Camera 

The  camera  contains  a  192  X  165  pixel  Texas  Instruments  TC21 1  CCD 
image  sensor.  Its  specifications  are  in  Appendix  A  [Ref.  19],  Figure  7  is  a 
photograph  of  the  sensor  centered  inside  the  camera.  The  other  components 
inside  the  camera  head  translate  the  clock  signals  from  the  interface  board  and 
amplify  the  CCD  output. 

Figure  8  shows  a  functional  block  diagram  of  a  TC21 1  CCD  image 
sensor  [Ref.  19].  Key  components  include  the  silicon  matrix  of  pixels  and  a  set 
of  gates  to  shift  the  charge.  The  image  area  gate  (IAG)  performs  a  parallel  shift 
of  all  rows  for  each  clock  pulse.  As  a  row  enters  the  output  serial  register  the 
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Figure  7.  Photograph  of  inside  the  CCD  camera  head  of  the  Spectra  Sources 
Lynxx  PC+  system. 

serial  register  gate  (SRG)  shifts  the  row  to  the  charge  detection  amplifier.  The 
TC21 1  includes  an  antiblooming  gate  (ABG)  that  should  not  be  used  since  it 
introduces  a  severe  nonlinear  photo  response. 


The  method  used  to  acquire  and  to  read-out  an  exposure  determines 
the  CCD  sensor’s  read-out  rate.  The  electrons  collect  in  potential  wells,  or  pixels 
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Figure  8.  Lynxx  PC+  image  sensor  CCD,  Texas  Instruments  TC211  block 
diagram. 

after  exposing  the  image  area  to  light.  The  charges  then  shift  down  to  the  output 
register  as  rows,  through  165  timing  cycles.  The  output  register  then  shifts  the 
charges  210  times  for  each  row.  This  cycle  of  210  shifts  includes  12  dark  pixels 
for  a  dark  reference  and  6  dummy  pixels  used  to  transfer  the  charges  out  of  the 
register.  Figure  9  shows  the  charge  transfer  process  [Ref.  19].  Another  factor 
to  consider  for  an  image  sensor  chip  is  the  noise.  The  noise  equivalent  signal 
for  the  TC211  is  150  electrons.  It  depends  on  (kTC)1*  /q,  where  C  is  the 
capacitance  of  the  read-out  charge  collector,  k  is  the  Boltzmann  constant,  T  is 
the  absolute  temperature  and  q  is  the  charge  of  an  electron. 

The  integration  time  needed  to  achieve  a  particular  signal-to-noise  ratio 
depends  on  both  the  telescope  optics  and  the  object’s  radiance.  It  will  be  a 
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Figure  9.  TC21 1  image  area  transfer  process,  uses  a  virtual  phase. 

constant  time  added  to  the  read-out  time  regardless  of  the  number  of  pixels  used 
or  transfer  times.  The  image  area  gate  (IAG)  and  serial  register  gate  (SRG) 
maximum  pulse  rates  determine  the  sensor  read-out  rate. 
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The  maximum  rate  for  the  IAG  parallel  row  transfer  gate  is  1.5  MHz. 
This  corresponds  to  0.667  microseconds  (jis)  per  row.  Clocking  the  IAG  shifts 
the  whole  array  down  by  one  row.  This  can  be  used  to  clear  the  CCD.  To  clear 
all  165  rows  requires  only  110  ^s.  Since  an  IAG  clock  pulse  clears  the  serial 
register,  the  IAG  should  not  be  clocked  any  faster  than  the  SRG  can  read-out  a 
row  of  pixels. 

The  maximum  rate  for  the  SRG  single  row  shift  gate  is  10  MHz.  This 
corresponds  to  0.1  ns  per  pixel.  Figure  8  shows  the  pixel  read-out  order.  Six 
dummy  pixels  will  appear  before  any  valid  data  and  the  twelve  trailing  dark  pixels 
can  be  ignored.  At  the  maximum  clock  rates,  the  ideal  minimum  time  required 
to  read  an  entire  row  of  pixels  is  210  pixels  X  0.1  \is  =  21  |is.  The  minimum  time 
needed  for  all  165  rows  is  165  rows  X  21  \l$  =  3465  \is. 

The  total  read-out  time  must  include  the  serial  shift  time  and  two 
parallel  array  shifts.  One  set  of  IAG  parallel  shifts  are  needed  to  move  each  row 
into  the  output  shift  register  and  a  second  set  must  clear  the  CCD  before  the  next 
exposure.  Since  the  output  register  cannot  be  read  during  a  row  transfer  (IAG 
clocking),  using  all  the  pixels,  the  ideal  minimum  read-out  time  is  3.685  ms. 
Placing  the  image  in  the  lower  left  portion  of  the  array  can  increase  the  frame  rate 
since  only  a  portion  of  the  array  needs  to  be  read-out.  After  shifting  the  rows 
with  information  to  the  output  register,  the  IAG  clocking  can  increase  to  its 
maximum  rate  to  clear  the  CCD  and  to  commence  another  integration.  The  six 
leading  dummy  pixels  will  always  add  some  time  to  each  serial  read-out  time,  but 
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the  twelve  trailing  dark  pixels  could  be  chopped  off.  As  an  example,  a  60  X  60 
image  would  take  2  X  60  rows  X  0.667  jiS  +  66  pixels  X  .1  jis  X  60  rows  for  the 
serial  read-out  for  a  total  of  476  jis.  There  are  66  pixels  per  row  because  of  the 
six  dummy  pixels. 

As  mentioned  earlier,  the  integration  or  exposure  time  adds  to  the  read¬ 
out  time  to  produce  the  total  measurement  time.  The  total  time  could  be  reduced 
to  only  the  integration  time  if  the  CCD  were  read-out  while  integrating  the  next 
set  of  data.  This  reduces  the  minimum  total  time  to  the  longer  of  these  two 
processes.  Frame  transfer  devices  as  shown  in  Figure  6  have  a  complete  extra 
set  of  pixels  to  hold  the  data  for  read-out.  The  interline  method  is  not  practical 
for  a  star  image  because  the  gaps  at  the  transfer  pixel  rows  reduce  the  exposed 
area  by  a  factor  of  two.  Each  pixel  would  lose  half  of  its  photons.  This  is 
unacceptable  for  this  application.  Masking  off  half  of  the  array  is  a  possible 
solution.  It  would  not  work  with  the  TC21 1  chip.  Clocking  the  IAG  parallel  row 
shift  register  would  produce  streaking  in  the  upper  half  exposed  pixels  while 
reading  the  masked  pixels. 

2.  Interface  Card 

The  Spectra  Source  computer  interface  card  contains  the  circuitry 
needed  to  transfer  information  from  the  CCD  array  to  the  computer  bus.  It  does 
an  eight  bit  parallel  transfer  to  an  IBM  compatible  computer.  Figure  10  is  a  block 
diagram  of  the  components  of  this  card.  The  CCD  sensor  was  discussed  in 
Section  1 .  The  sample-and-hold  (S&H)  amplifier  chip  converts  video  voltage  from 
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Figure  10.  Spectra  Sources  Lynxx  PC+  interface  card  block  diagram. 

Section  1.  The  sample-and-hold  (S&H)  amplifier  chip  converts  video  voltage 
from  one  pixel  into  a  voltage  pulse.  It  is  a  NEC  SE/NE5537  [Ref.  20].  The 
minimum  acquisition  time  is  less  than  4  jis.  It  is  controlled  by  pulses  from  the 
chip  below  it,  a  74LS123P  [Ref.  21],  which  is  a  dual  retriggerable  one-shot.  Its 
pulse  switching  characteristics  are  in  the  nanosecond  (ns)  range.  A  1200  pF 
capacitor  and  1 2  k  resistor  determine  the  present  pulse  length  of  51 80  ns.  The 
control  chip  characteristics  did  not  limit  the  speed  of  the  S&H  chip. 

The  analog  digital  converter  (ADC)  chip  converts  the  voltage  to  a 
positive  12  bit  integer,  which  has  a  maximum  value  of  4095.  This  system  has  a 
Maxim  MX7572  ADC  chip  [Ref.  22].  Its  minimum  conversion  time  is  5  *is.  It,  as 
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and-hold  and  digital  conversion  time  is  9  ^s  for  each  pixel.  For  a  60  X  60  pixel 
image,  the  total  S&H  and  ADC  time  is  32.4  ms. 

3.  Computers 

Two  different  computers  were  used  to  investigate  the  dependence  of 
the  frame  rate  on  computer  characteristics.  The  first  was  a  Compaq  386,  20 
megahertz  machine.  It  had  a  16  megabyte  RAM  and  a  60  megabyte  hard  disk. 
The  second  was  a  Dell  486,  33  megahertz  machine.  It  had  a  16  megabyte  RAM 
with  a  350  megabyte  hard  disk. 

B.  SOFTWARE 

The  Lynxx  system  came  with  software  programs  designed  for  amateur 
astronomical  observations  [Ref.  17].  This  thesis  used  a  separate  set  of  data 
acquisition  program  modules  written  for  NPS  by  Spectra  Sources  [Ref.  18].  It 
was  necessary  to  modify  these  modules  during  this  thesis  research.  The 
software  was  written  in  the  C  language.  Microsoft  Quick  C  2.5,  Borland  Turbo 
C++  1.0  and  Borland  C++  2.0  compilers  produced  executable  code  that  were 
tested  to  compare  their  relative  speeds  while  collecting  and  processing  data. 
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IV  RESULTS 


The  first  task  involved  writing  the  modulation  transfer  function,  centroid  jitter 
and  display  programs.  After  it  become  clear  that  the  Lynxx  PC  +  CCD  system 
was  too  slow  for  effective  centroid  measurements,  the  second  task  developed 
programs  to  analyze  the  Lynxx  image  acquisitions  times.  The  data  acquisition 
functions  were  timed  using  different  computers  and  code  from  different 
compilers.  The  manufacturer’s  maximum  speed  specifications  for  both  the  CCD 
camera  and  the  interface  card  were  compared  to  the  actual  measured  times. 
After  analyzing  these  results,  replacements  were  recommended. 

A.  PROGRAMS 

A  modulation  transfer  function  program  and  a  jitter  program  were  written  to 
compute  the  atmospheric  coherence  length,  r0.  An  artificial  star  program  was 
developed  to  test  the  MTF  program.  The  following  sections  describe  these 
programs. 

1.  Modulation  Transfer  Function 

This  program  calculated  the  coherence  length,  r0,  from  a  stellar  image 
intensity  distribution.  It  subtracted  the  background  from  the  star  image, 
calculated  the  line  spread  function,  found  the  centroid  and  calculated  the  fast 
Fourier  transform.  It  then  calculated  r0  using  equation  (9).  The  calculations  were 
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done  in  the  x  and  y  dimensions  separately,  giving  a  quality  check,  since  actual 
image  turbulence  is  nearly  always  isotropic.  With  the  image  centroid  data  the 
telescope  could  operate  in  an  auto  track  mode,  a  feature  that  comes  with  the 
Lynxx  software  package.  The  program  plotted  the  modulation  transfer  function 
as  a  visual  check.  The  program  separates  into  modules  for  ease  of  future 
revisions.  This  program  is  in  Appendix  B. 

2.  Artificial  Star  Program 

This  program  created  a  perfect  exponential  star  image  for  the  MTF 
program.  This  provided  a  convenient  test  of  the  MTF  program,  since  the  FFT  of 
an  exponential  is  also  an  exponential.  This  program  is  in  Appendix  C. 

3.  Jitter 

This  program  calculated  rQ  from  the  centroid  jitter.  It  first  removed  the 
background  by  subtraction.  After  this,  it  computed  the  image  centroid  jitter  using 
the  standard  deviation  of  the  image,  according  to  equation  (11).  This  program 
also  calculated  the  line  spread  functions  and  plotted  them  on  the  screen  so  that 
proper  telescope  tracking  was  verified.  Appendix  D  contains  this  jitter  program. 

B.  TIME  MEASUREMENTS 

Programs  were  written  or  modified  to  compare  frame  rates  attainable  for 
different  subframe  sizes.  The  computer  clock  measured  the  time  differences. 
Since  the  MS-DOS  system  clock  had  a  55  ms  resolution,  many  of  the  times  were 
done  for  1 00  cycles  of  the  function  of  interest.  All  the  times  used  were  at  least 
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one  second  long,  resulting  in  an  accuracy  within  1%.  An  array  of  5  or  more  trials 
ensured  consistency.  Appendix  E  is  an  example  of  the  code  used  for  measuring 
the  times. 

1.  Computer 

We  found  that  the  type  and  speed  of  the  computer  did  not  contribute 
significantly  to  the  camera  frame  rate.  Measurement  time  included  the  clearing 
of  the  CCD,  the  integration,  and  digitization  times.  With  an  image  size  of  60  X 
60,  there  was  less  than  Vs  of  a  frame  per  second  difference  between  386  and 
486-based  computers.  Figure  1 1  shows  a  comparison  of  frame  rates  verses 
subframe  size  for  the  two  machines.  Including  calculations  for  the  line  spread 
functions  and  backgrounds,  the  frame  rate  improvement  was  about  3Vs  frames 
per  second  for  a  60  X  60  image  with  a  486  computer  over  that  achievable  with 
a  386  machine.  This  is  shown  in  Figure  12.  The  mathematical  calculations  were 
timed  to  determine  their  contributions  to  the  total  frame  rate.  These  calculations 
could  have  been  done  separately  if  they  had  slowed  the  frame  rate,  but  they 
were  actually  performed  in  real  time. 

2.  Compiler 

It  was  found  that  the  choice  of  compiler  did  not  significantly  influence 
the  speed  of  the  measurements.  Figure  13  shows  there  was  no  difference  in  the 
frame  rate  for  the  code  produced  by  three  different  compilers  to  acquire  data. 
Figure  14  shows  that  after  adding  the  mathematical  calculations  to  the  acquisition 
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COMPARISON  OF  PROCESSORS  WITHOUT  CALCULATIONS 
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Figure  11.  Spectra  Sources  Lynxx  PC+  frame  rates  using  IBM  PC 
386  and  486  computers.  Times  include  clear,  integration  and 
digitization,  using  code  compiled  under  Turbo  C++. 
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time,  the  code  produced  by  Quick  C  was  about  %  of  a  frame/second  slower  than 
the  code  from  Borland  C++  or  Turbo  C  +  +  for  a  60  X  60  image  size. 

3.  Functions 


The  times  for  each  measurement  were  divided  into  four  steps.  Table 
1  shows  these  steps. 

TABLE  1 


Step 

Image  size 

Time 

Clear  CCD 

Must  clear  entire  array 

constant  6.45  ms 

Integration 

Independent  of  s'ze 

used  5.0  ms 

Digitization 

60  X  60 

86.1  ms 

Calculations 

60  X  60 

40.2  ms 

Total 

60  X  60 

137.75  ms 

The  digitization  and  calculation  times  depended  on  the  size  of  the  image,  while 
the  clear  and  integration  times  were  constant  for  each  measurement.  This  is 
shown  in  Figure  15.  Clear  and  digitization  times  depended  on  the  image  sensor 
and  interface  card. 

C.  HARDWARE  RATE  COMPARISONS 

This  section  compared  the  maximum  design  speeds  of  the 
chips  in  both  the  camera  head  and  interface  card  to  those  actually  measured. 
A  general  time  reference  was  200  frames/second  which  corresponded  to  5  ms 
per  frame. 
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1.  CCD  Camera 


The  measured  time  for  the  clear  function  was  6.45  ms.  The  present 
sensor,  TC211,  is  capable  of  clearing  the  entire  array  in  110  jis.  All  rows  need 
to  be  cleared  even  when  reading  a  sub  array  such  as  a  60  X  60  image.  This 
function  was  much  slower  than  the  sensor  capability.  It  is  not  known  why  the 
clear  times  are  so  slow. 

The  present  sensor  could  read  a  60  X  60  image  in  436  jis  at  maximum 
design  speeds.  The  measured  digitization  time,  which  included  the  read-out 
time,  sample-and-hold  (S&H)  and  the  digital  conversion  (ADC)  times  was  86.1 
ms.  Since  this  was  much  larger  than  436  ns,  the  acquisition  time  depended 
primarily  on  the  S&H  and  ADC  processes. 

2.  Interface  Card 

The  interface  card  controls  the  S&H  and  ADC  times.  Their  maximum 
speed  specifications  would  result  in  processing  a  60  X  60  image  in  32.4  ms.  This 
is  a  significant  portion  of  the  observed  86.1  ms  digitization  time.  We  therefore 
concluded  that  these  two  chips  were  the  primary  contributors  to  the  frame  rate 
and  were  not  operating  at  their  maximum  speed  specifications.  The  sensor  itself 
can  hypothetically  operate  at  significantly  higher  speeds. 
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D.  RECOMMENDATIONS 


1.  CCD  Camera 

To  increase  the  frame  rate  requires  an  image  sensor  that  uses  a  frame 
transfer  method  with  separate  image  and  storage  arrays.  The  Texas  Instruments 
TC277,  a  735  X  580  pixel  CCD  image  sensor  is  a  proposed  replacement  sensor 
{Ref.  19],  Figure  16  shows  its  set-up.  The  image  area  is  699  X  288  which  is  more 
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Figure  16.  Recommended  image  sensor  CCD,  Texas  Instruments  TC277  block 
diagram. 

than  big  enough  for  our  applications.  The  image  sensing  area  and  image 
storage  area  have  different  gates.  The  image  area  gate  (IAG)  shifts  the  rows  from 
the  sensing  area  to  the  storage  area  one  row  per  clock  pulse.  The  storage  area 
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gate  (SAG)  shifts  the  rows  down  to  the  serial  register  at  one  row  per  clock  pulse. 
The  three  serial  register  gates  (SRG)  shift  out  the  pixels  in  each  row. 

The  storage  area  gate  (SAG)  maximum  clock  rate  is  3.34  MHz,  which 
corresponds  to  0.299  ns  per  row.  The  serial  registers  gates  (SRG)  maximum 
clock  rates  are  4.46  MHz,  which  corresponds  to  0.224  ^s  per  pixel.  For  a  60  X 
60  image,  noting  that  there  are  eleven  dummy  pixels  per  row,  the  read-out  time 
is  60  rows  X  0.299  ^s  +  71  pixels  X  0.224  *is  X  60  rows  for  a  total  of  972  iiS. 
This  is  done  in  parallel  with  the  integration.  Since  this  is  much  less  than  any 
expected  integration  time,  the  integration  time  and  transfer  time  of  data  from  the 
sensing  area  to  the  storage  area  become  the  limiting  factors  for  frame  rates. 

The  image  area  gate  (IAG)  maximum  clock  rate  is  3.34  MHz,  which 
corresponds  to  0.299  ns  per  row.  The  clear  time  would  then  be  288  rows  X 
0.299  ^s  or  86  \is.  This  clear  time  plus  the  number  of  rows  X  0.299  fis  plus  the 
integration  time  will  be  the  measurement  time  for  each  set  of  data.  For  an  image 
with  60  rows,  the  measurement  time  would  be  only  86  +  17.9  =  104  jis,  plus  the 
integration  time.  For  the  TC21 1 ,  the  minimum  clear  and  60  X  60  image  read-out 
time  is  110  +  436  =  546  jis.  The  TC277  is  over  5  times  faster. 

It  should  also  be  noted  that  the  number  of  pixels  in  each  row  of  the 
TC277  can  be  the  maximum  available  (735)  without  affecting  the  speed,  due  to 
the  parallel  read-out  and  integration  times.  The  noise  equivalent  signal  is  only 
25  electrons,  a  6-fold  improvement  from  the  TC211’s  150  electrons.  Replacing 
the  image  sensor  would  require  modifications  to  the  interface  board. 
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Using  such  an  improved  sensor,  the  bulk  of  the  data  could  be 
transferred  while  collecting  the  next  data  set.  What  remains  is  to  get  the 
sampling  and  digital  conversion  times  fast  enough  to  be  completed  during  the 
integration  time. 

2.  Interface  Card 

The  first  step  in  the  digital  conversion  is  the  sample-and-hold  (S&H). 
A  proposed  replacement  chip  for  this  is  Analog  Devices  AD781  [Ref.  23],  It 
provides  a  700  ns  conversion  time.  This  is  over  5  times  faster  than  the  present 
SE/NE5537  S&H  chip’s  4  ns.  For  a  60  X  60  image,  the  AD781  S&H  time  reduces 
to  2.52  ms  from  14.4  ms  for  the  present  chip.  The  present  control  chip, 
DM74LS123P,  for  the  S&H  chip  does  not  need  replacement.  Changing  its 
capacitor  to  55  pF  and  resistor  to  25  k  would  adjust  its  pulse  width  to  700  ns 
[Ref.  21].  This  is  the  only  revision  necessary  to  replace  the  S&H  chip. 

A  proposed  replacement  chip  for  the  ADC  is  a  Maxim  Maxi  62  [Ref.  22]. 
This  chip  is  almost  identical  to  the  present  one,  except  that  its  conversion  time 
is  3  *is  instead  of  5  jxs.  For  a  60  X  60  image,  the  improved  ADC  time  would  be 
10.8  ms.  The  Maxi  62  could  replace  the  Mx7552  with  little  revision  to  the  rest  of 
the  circuit. 

There  are  some  faster  S&H  and  ADC  chips  available  if  the  board  were 
rebuilt.  Datel  has  a  S&H  chip  with  a  time  of  only  25  ns  and  a  ADC  chip  with  a 
conversion  time  of  350  ns  [Ref.  24].  A  better  option  is  to  have  the  S&H  and  ADC 
processes  combined  in  one  chip.  Of  those  available,  the  Datel  ADS1 18  [Ref.  24] 
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has  a  throughput  rate  of  5.0  Mhz,  which  corresponds  to  only  200  ns  per  pixel. 
The  DM74LS123P  control  chip  will  still  work  by  using  a  60  pF  capacitor  with  a  5 
k  resistor  [Ref.  21].  Besides  providing  fast  timing,  combining  these  two  chips 
would  also  simplify  the  surrounding  circuit.  For  a  60  X  60  image,  the  S&H  and 
ADC  time  would  be  720  jis  for  the  ADS1 18. 

E.  SUMMARY 

Table  2  summarizes  all  of  the  relevant  times.  The  components  are  replaced 
going  from  left  to  right  across  the  table.  They  are  in  order  of  the  complexity  of 
replacement.  Replacing  the  CCD  sensor  would  be  the  significant  step  requiring 
major  modifications  to  the  interface  board.  There  are  also  parallel  times  after  the 
CCD  sensor  is  replaced.  The  processes  of  clearing,  integration  and  image  area 
gate  (IAG)  transfer  occur  parallel  to  the  read-out  and  digital  conversion.  The 
slowest  of  these  times  will  limit  the  frame  rate.  The  S&H  and  ADC  times  are  the 
limiting  factors.  For  this  reason  and  the  necessity  to  modify  the  board  after 
replacing  the  CCD  sensor,  it  is  best  to  build  a  new  system. 
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'The  lowest  of  these  two  rates  will  limit  the  system.  These  two  rates  reflect  the  parallel  times  introduced  by  the 
frame  transfer  CCD. 


V  CONCLUSIONS  AND  RECOMMENDATIONS 


A  CCD  camera  can  measure  atmospheric  turbulence  by  measuring  the 

image  spread  of  a  point  source  and  by  measuring  the  centroid  jitter  induced  by 

turbulence.  The  Spectra  Sources  Lynxx  PC+  system  evaluated  in  this  thesis  can 

measure  the  coherence  length  using  the  point  source  image  spread  technique, 

but  it  is  much  too  slow  for  jitter  measurements.  Centroid  motion  and  Greenwood 

frequency  measurements  need  a  sample  rate  of  200  Hz  or  more  to  avoid 

undersampling  the  atmospheric  dynamics.  Achievement  of  these  rates  requires 

a  frame  transfer  CCD,  so  that  parallel  image  exposure  and  frame  read-out  are 

possible.  Data  processing  can  be  performed  during  an  exposure  if  the  time 

needed  is  sufficiently  short,  or  after  digitizing  a  series  of  exposures.  The  Lynxx 

PC+  system,  which  provides  12  frames/second,  is  much  too  slow  for  this.  Since 

both  the  CCD  sensor  and  digital  conversion  components  need  replacement,  it 

is  best  to  design  and  build  a  new  system.  The  recommended  CCD  image  sensor 

is  a  Texas  Instruments  TC277  [Ref.  19].  The  recommended  replacement  for  the 

sample-and-hold,  and  analog  digital  converter  chips  is  a  combined  sampling 

analog-to-digital  converter  chip,  Datel  ADS118  [Ref.  24].  The  atmospheric 

« 

coherence  length  programs  developed  in  this  thesis  should  work  with  a  new 
system. 
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APPENDIX  A  TEXAS  INSTRUMENTS  TC211  CCD  IMAGE  SENSOR 


TC211 

192  x  165  PIXEL  CCD  MAGE  SENSOR 

_  _  Ptt/t.  JMjUMjy  1*90 


•  Full F  ram*  Operation 

•  Anilblooming  Capability 

•  Single  Phase  Clocking  for  Horizontal  and 
Vertical  Transfers 

•  Fast  Clear  Capability 

•  Dynamic  Range . . .  60  dB  Typical 

•  High  Blue  Response 

•  High  Photoresponse  Uniformity 

•  Solid  Stale  Reliability  With  No  Image 
Burn  In,  Residual  Imaging,  Image 
Distortion,  Image  I  ag,  or  Microphonics 

•  6  Pin  Dual  In  line  Ceramic  Package 

•  Square  Image  Area: 

-  2640  |im  by  2640  itm 

-  192  Pixels  (II)  by  165  Pixels  (V) 

-  Each  Pixel  13.75  pm  (11)  by  16  pm  (V) 

description 

The  1C2I1  is  a  fuN  Irame  charge  coupled  device  (COO)  Image  sensor  dosigned  specifically  lor  industrial 
applications  requiring  iiiggedness  and  small  sire  1  he  Image  sensing  area  is  configured  krto  IBS  horizontal  lines 
each  containing  192  pivots  Iwolvo  additional  pixels  are  providud  el  tlie  end  ol  each  line  to  establish  a  dark 
relerence  and  line  clamp  Tlie  aniMootning  feature  Is  activated  by  supplying  clock  pulses  lo  the  anlibtoomkvg 
gale,  an  Integral  pari  ol  each  Image  sensing  element  The  charge  is  converted  lo  signal  voltage  at  4  pV  per 
electron  by  a  high  performance  slnrcture  with  built  In  automatic  rosol  and  a  voltage  reference  generator.  The 
signal  is  lurther  bulleted  by  a  low  noise  two  stage  source  tottowor  amplifier  lo  provide  high  otrlpnl  drive 
capabthty 

The  TC211  Is  supplied  In  a  B  pin  dual  In  line  ceramic  package  approximately  7,5  mm  (0.3  In  )  square  and  la 
characterized  lor  operation  koni  -  10‘C  to  45  “C  Tire  glass  wktdowcan  be  cleaned  using  any  standard  method 
for  cleaning  optical  assemblies  or  by  wiping  lire  sntluco  with  a  cotton  swab  soaked  in  alcohol. 
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TC211 

192  x  165-PIXEL  CCD  IMAGE  SENSOR 
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functional  block  diagram 


pin  functional  description 
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functional  description 

The  Image  sensing  atea  conslsls  ol  185  liorizontal  Image  Hues  endt  containing  192  ptiolosenslttve  elements 
(pixels)  E  acli  pixel  Is  1 375  pm  (horizontal)  by  1 8  00|im  (vertical)  As  Hghl  enters  lire  silicon  In  Hie  Image  sensing 
area,  free  electrons  are  generated  and  collected  In  potential  weds  (see  rignre  1)  During  (its  time.  Hie 
antiblooming  gate  Is  activated  by  applying  a  burst  ol  pulses  every  horizontal  blanking  Interval  Tills  prevents 
blooming  caused  by  the  spilling  ol  charge  Irom  overexposed  elements  kilo  netgldjorlng  elements.  The 
antiblooming  gate  Is  typically  held  at  a  mid  level  voltage  di  sing  readout.  The  quantity  ol  charge  coltected  In  each 
pixel  Is  a  linear  function  of  tie  incident  dgtrt  end  lie  exposure  lime.  After  exposure  and  under  dark  condHIons, 
Hie  charge  packets  are  transferred  from  Hie  Image  area  to  Hie  serial  register  at  tie  rate  ol  one  fctiage  Rue  per 
each  dock  pulse  applied  to  Hie  Image  area  gale  Once  an  Image  due  lias  been  transferred  Into  Hie  serial  register. 
Hie  serial  register  gale  can  be  clocked  unld  ad  ol  Hie  dtarge  packets  are  moved  out  ol  the  serial  register  to  Hie 
diarge  detection  node  at  Hie  amplifier  Input 

There  are  12  dark  pixels  to  the  right  oltlie  192  Imago  pixels  on  eadi  Imago  Hue  These  dark  pixels  are  shielded 
Irom  Incident  dgfit  and  the  signal  derived  Irom  fiiem  can  be  used  lo  generate  a  dark  reference  lor  restoration 
ol  ttie  video  black  level  on  Hie  next  Image  due 
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TC211 

192  x  165  PIXEL  CCD  IMAGE  SENSOR 
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Each  clock  pulse  applied  to  Die  image  atoa  pale  causes  an  automatic  Iasi  cleat  ol  Die  t#2  Image  pixels  and  12 
dark  pixels  ol  the  serial  register  before  Ihe  next  image  Sue  Is  Itansfoued  into  Die  serial  regisler  (Note  foal  foe 
six  dummy  pixels  al  llw  front  ol  Die  serial  regisler,  which  are  itsod  lo  transport  charge  packets  from  foe  serial 
regisler  to  Die  ampWter  input,  are  not  cleared  by  Die  Image  aroa  gale  dock )  The 'automatic  test  dear  feature 
can  be  used  to  friWaNze  lie  Image  area  by  transferring  al  16S  Image  Ines  lo  Die  serial  regisler  gate  under  dark 
conditions  witliout  docking  foe  serial  rogisler  gate. 
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Figure  1.  Charge  Accumulation  and  Transfer  Process 
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TC211 

192  x  165  PIXEL  CCD  IMAGE  SENSOR 

03?/r,  JANUARY  1990  _  _ 
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Figure  2.  liming  Diagram,  Nonlnietlace  Mode 


absolute  maximum  ratings  over  operating  tree-air  temperature  range  (unless  otherwise  noted) 


Supply  voltage  range  (ADO)  (see  Nolo  1)  .  0  V  lo  15  V 

Clock  voltage  range  lor  IAG.  SflG.  ABG  terminals  .  -15Vlo5V 

Operating  tree  air  temperature  range  .  -  30'C  lo  85*C 

Storage  temperature  range  .  -30~CtoBS*C 

I  end  temperate  e  t.8  mtn  (t/lfl  Inch)  Ironi  case  lor  10  seconds  .  260"C 
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192  x  165  PIXEL  CCD  IMAGE  SENSOR 

_  _ ran. January  two 
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TC211 

192  x  165-PIXEL  CCD  IMAGE  SENSOR 


P3271.  JANUARY  IMP _ * 


electrical  characteristics  over  recommended  operating  range  of  supply  voltage,  TA  =  -  10"C  to  40X 


PARAMETER 

mm  mi  max 

UNIT 

Dynamic  range  (see  Mole  2) 

AnttAoomlng  dtenUod  («e»  Note  3) 

60 

dB 

AntMuorrfng  enabled 

57 

Charge  conversion  (actor 

mmamnmm 

|VM 

Charge  tr antler  eHictency  (tea  Note  4) 

^H^SZZZHHni 

Signal  response  delay  Brno, » (taa  Note  8  and  Figure  5) 

mam 

Qwnmi  |M*  Mote  t) 

SB 

OiRpul  tnManca 

iMHMMClBiRKZJ 

mm 

Noise  vofage 

1/1  noise  (S  kHi) 

J70 

nW/TE 

Random  noise.  I  -  100  Wt* 

|  Nolee  equivalent  tignal 

huh 

Refaction  ratio  at  7  16  MM* 

From  A 00  to  output  (toe  Note  7) 

fi  KA 

n 

From  SRO  lo  output  (toe  Note  8) 

37 

|  Stpply  currant.  »oo 

5  10 

Capacitance 

Image  area  gate 

«M0 

H 

Sofia!  register  gale 

25 

rssBaammmmmm 

I  AJI  typical  valuei  eta  el  T*  -  25'C 


MOT  f  S:  2  Dynamic  r ange  It  - 20  timet  Hie  logmHfim  ol  Hie  mom  twite  t«gnnf  dh/«fnr1  !>y  Hie  tnfureHon  outprrf  tJgnnl 

3  For  H4s  leal,  H»e  mill  dooming  gate  mutt  lm  tiinsnd  at  Hie  iMtormcdtafo  loval 

4  Charge  liemlet  efficiency  It  one  ml™  is  H*e  rl*ar go  lot®  per  transfer  lo  H*o  output  register  Tim  teat  It  portormod  In  the  dwfc  uelug  an 
electric  nl  lo|Mit  signal 

$  Signal  response  delay  time  It  Vie  Mine  between  Hie  Infing  rxlga  ot  Hie  SIX*  crock  pohe  end  Itie  output  signal  veSrt  tlola 
6  Gemma  (y)  fo  lie  vntue  ot  fie  evpnoont  In  the  equation  below  lot  two  {minis  on  lie  Street  portion  ot  H«e  transfer  function  curve  (Nt 
vnfcie  represents  point®  neat  tab  it  af  on) 

(Exposure  (2)\r  _  / Qutprd  signal  (?) 

Exposure  (1)  J  ^Output  signal  ft) 

7  AD8  toJocUon  tnlio  It  -  20  timet  H*e  tognHImi  ol  Hkht  nmpttu'lo  nt  Hie  out|Nit  drvitlndby  H»o  nr  arnfAtudo  at  A 00 
8  SRO  rafncfon  ratio  It  -?0  ffmot  H*e  tngnrifim  ol  tl*o  ec  anyJWiKlo  nl  ll*e  output  rfvirfnrf  tiy  H»e  ac  ampMode  at  SRO 


,  IfeXAS  ^ 
Insihuments 

2-  li  rost  omct  eo«t«5w*0NiM.  itxAinm 
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TC211 

192  x  165-PIXEL  CCD  IMAGE  SENSOR 

_ _ _ P»7|.  JANUARY  WOO 


optical  characteristics!  =  25  C  (unless  otherwise  noted) 


OGJZ3 

EH 

StMisNivay  (see  Note  0) 

C3I3I3H 

260 

mVAa 

33 

Saba  abort  signal  (see  Note  11) 

AnbUoomlng  (ftaaNod 

4oo  eoo 

AnIWooming  enaUod 

pEZ^KZSH 

Blooming  overload  lato  (toe  Note  12) 

Sfcuhe 

plHHI 

Shuttorodfegld 

nUMESAl 

Output  signal  nonuntformity  (1/2  saturation)  (see  Note  16) 

10%  20% 

Image  area  wal  capacity 

ISOxtO91 

Dark  current 

|fA-*l*C 

0027 

Dark  signal  (see  Note  14) 

D»k  signal  nonunifcfmity  to.  sntss  laid  (so*  Noia  IS) 

imiini  ■  n  i 

mV 

ModuLilron  Ira  fib  let  ktocbon 

Horizontal 

60% 

70% 

NOTES:  0  Sonsalivrly  ie  measured  el  an  integialion  lime  of  16  667  me  and  a  source  lomporakite  of  2656  K  A  CM  S00  Rlor  it  used 
10  V||  it  ftw  output  voSag«  fiat  represents  Ihe  Rtf  oslioM  of  of*w  abort  of  ariMJnumeiy  V|j  •.  1/2  % ake <4ion  signal 
If  Saturation  it  fit  condition  In  which  (orthoi  inciausa  In  tiposiwo  duos  nut  load  to  bettor  bKfoaso  btoutfNj  wqimI. 

12  8tonrning  ovtiiload  ratio  it  tttu  toko  of  Hoomirtg  eiposoro  lo  tutur.ilion  o*pnsote 

13  Output signal nonontformity  btt>e » alio olRtonta »lmi impisuf  lo  |«1aot tklor our  o  in  output  signal  to  N>o  mean  o^mtstgnd  lot  eipoatee 
O'V  islod  lo  give  1/2  Itiu  totm abort  otitpul  tign.il 

14  Oatk  signal  level  »  moamrel  bom  tie  dummy  (null 

15  Dark  signal  noniNMlotmtly  is  Hm>  manmum  |x»d  lo  piavl  ikluience  in  a  d.*k  corvfcfcon 


PARAMETER  MEASUREMENT  INFORMATION 


Figure  3.  Typical  Clock  Waveform  for  IAG  and  ABG 

SIEWRAIE  BETWEEN  T 0%  and  00%  -  70  to  120  W|.a  •,  -  ISO  na.  I|  -  SO  ns 
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TC211 

192  x  165  PIXEL  CCD  IMAGE  SENSOR 


Ptt/I,  JANUARY  I  WO 


PARAMETER  MEASUREMENT  INF' 


Figure  «.  Typical  Clock  Waveform 


SAMPIf 

AND 

MOlD_ 


Figure  5.  SRO  end  CCD  Output  Wai 


# 


2 20 
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TC211 

192  x  165-PIXEL  CCD  IMAGE  SENSOR 


mj  IXC71,  JANUARY  WM 

TYPICAL  CHARACTERISTICS 


HORIZONTAL  MODULATION 
TRANSFER  FUNCTION 
(BARS  PARAl  LEL  TO  SERIAL  REGISTER) 


vertical  modulation 
TRANSFER  FUNCTION 

(BARS  PERPENDICULAR  TO  SERIAL  REGISTER) 


NormaNitd  Spatial  Fraquancy 
I _ I _ I _ I _ I _ I 

o  n  us  is  s  iso  ns 

Spatial  Praquancy  -  Cyclaafaim 


Norm aftxad  Spatial  praquancy 

I _ I _ I _ I _ I _ I 

0  M  l«t  Ilf  Ml  SOS 

Spatial  Praquancy  -  CydaaMww 


Figure  fl 


Figure  1 


NOISE  SPECTRUM  OF  OUTPUT  AMPLIFIER 


CCD  SPECTRAL  RESPOMSIVITY 


lock! ml  Wrvrlmglli  -  nm 


Figure  • 


Figure* 
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TC211 

192  x  165  PIXEL  CCD  IMAGE  SENSOR 


03271,  JANUARY  I  WO 

TYPICAL  APPLICATION  DATA 


Figure  10.  Typical  Application  Circuit  Diagram 


,  TfeXAS  ^ 
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TC211 

192  x  165  PIXEL  CCD  IMAGE  SENSOR 

_ _ DM>I.J*HLI*WT  IMP 


MECHANICAL  DATA 

* 

T  lie  package  lor  the  TC2II  consists  ol  a  ceramic  base,  glass  window,  and  a  8  lead  frame.  The  glass  window 
is  sealed  to  (lie  package  by  an  epoxy  adhesive.  The  package  luads  are  configured  in  a  dual  in  Nne  organization 
and  lit  into  moimting  holes  with  2,54  mm  (0  I  Inch)  centerto  center  spacirtgs 


NOTES  A  Dhiiuomoiis  ate  in iiHtmtoluta  •iKjpa'ottMtttlK.oHy  hi  hh  Ih>*  Srnylo  Umtenwotw  «m»  nomniiri 
B  Tito  conic*  of  lt»o  |>a<  fuxju  aiv.!  t»»o  coulot  of  Uto  imago  ixoa  bio  im4  c  flint  ttlutif 

C  T  tio  cA^I.hko  lium  It  to  loft  ol  Ihe  gl.tbb  In  llto  itrt.tyo  sottsc**  ymf.tco  ts  l/|»*  nNy  1ntinf0  040iiK.il)  TliugluMMty|)icuVyO  Q20feidi1l4cfc 
b»mI  l»ns  on  hkIoi  of  tuft  action  of  1  52 
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APPENDIX  B  MTF  PROGRAM 


/*MTF.c  find  ro  and  f  from  the  Isf  using  FFT  to  the  MTF*******/ 

#define  TITLE  "Measure  ro,  plot  the  MTF" 


#define  AUTHOR 

/***********QyjQ|( 


"By  W.  J.  Rail" 

Q ******** *★*★★********★★★★★*★***★*★/ 


/*#include  <alloc.h>  */ 
#include  <lynxx.h> 


#define  True  1 
#define  False  0 
#define  false  0 


#include  <stdio.h> 

#include  <stdlib.h> 

#include  <graph.h> 

#include  <math.h> 

#include  <sys\types.h> 

/*int  frame_ptr[FRAME_BYTES/2];*/ 
int  size =60,  xstart=0,  ystart=0; 
float  lsfx[165],lsfmax=0.0,lsfy[192]; 

int  sign=1,  m=8,  mm=2,  mm2;/*mm  is  #  of  points  in  FFT,  mm2=mm/2***/ 
fl oat  mtfx[ 1 024] ,  mtf max = 0.0,  mtfy [1024],  re [ 1 024] , im  [ 1 024] , px[51 2] , py [51 2] ; 
float  rocal  =  1.0,rox,roy; 

void  init(); 

void  subexpose(long  tm,int  shutter); 
void  expose(long  tm); 
void  long_delay{long  tm)  ; 
void  calc Jsf (unsigned  int  far  *sub_frame); 
void  background(); 

void  fft(int  m,  int  sign,  float  re[],  float  im[]); 
void  calc_ro(); 
void  init_mtf(); 
void  plot_mtf(); 

/*to  see  Isf’s,  use  these  two  functions  instead  of  mtf  ones*/ 
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void  init_graph(); 
void  plot_graph(); 


/* - main  routine - 

void  main() 

{ 

int  i; 

for  (i=0;  i<m-1;  i++)  mm  =  2*mm; 
initJynxxO; 
init(); 

cooler_off(); 

_setvideomode(  _TEXTC80); 

> 

/* - Initialize  input  info - 

void  init() 

{ 

int  cool  =  0,shutter=0; 
long  tm; 

printfC'Expose  time  in  ms.:"); 
scanf("%ld",&tm); 
printf("Cooler  0=OFF  1=ON:“); 
scanf("%d,,,&cool); 
if  (cool) 

cooleronO; 

printf("Use  shutter  ?  (1  =Yes/0=No)"); 
scanf("%d",  &shutter); 
printf("Subframe  size?  "); 
scanf("%d",  &size); 

/*  printf("\nHit  any  key  to  quit:");*/ 

subexpose(tm, shutter); 

} 

/* - subExpose  CCD - 

void  subexpose(long  tm.int  shutter) 

{ 

unsigned  int  *sub_frame; 
int  sub_frame_bytes; 

float  pxO,pyO; 
int  i.j.x.y; 

for  (i=0;  i<m-1;  i++)  mm  =  2*mm; 


mm2 =(int)  ((float)  (mm)/2.0); 


sub_frame_bytes  =  size*size  *  sizeof(unsigned  int); 
sub_frame  =  (unsigned  int*)malloc(sub_frame_bytes); 

if  (!  shutter)  open_shutter(); 

while  (  IkbhitO  ) 

{ 

if  (shutter) 

{ 

expose  (tm); 

} 

else 

{ 

clrccd();  /*  if  not  using  shutter  *****************/ 
long_delay(tm); 

} 

xparam  =  (unsigned  char)  xstart; 
yparam  =  (unsigned  char)  ystart; 
size_param  =  (unsigned  char)  size; 
ptrparam  =  subframe; 
digitize_sub_frame  () ; 
calclsf  (subframe) ; 
background  (); 

for  (x = xstart + size- 1;  x>=xstart;  x--)  re[x]  =  lsfx[x]; 
fft(m,  sign,  re,  im); 

for  (x=0;  xcmm;  x++)  mtfx[x]  =  re[x]; 

/*find  the  x  power  Spectrum***************/ 
pxO  =  re[0]*re[0]  +  im[0]*im[0]; 
for  (j=0;  j<mm2;  j++) 

{ 

if(mtfx[j]  !=  0.0  )  if(px0  !=  0.0) 

px[j]=sqrt((re[j]*reO]+im[j]*imO])/pxO)/mtfx[j]; 

} 

/♦reinitialize  the  FFT  arguements  before  finding  the  Y  FFT*/ 
for  (i=0;  i<  1 024;  i++) 

{ 

re[i]  =  0.0; 
im[i]  =  0.0; 

} 
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for  (y =ystart+ size-1 ;  y>=ystart;  y--)  re[y]  =  lsfy[y]; 
fft(m,  sign,  re,  im); 

for  (y=0;  y<mm;  y++)  mtfy[y]  =  re[y]; 

/*find  the  y  power  spectrum***************/ 
pyO  =  re[0j*re[0]  +  im[0]*im[0]; 
for  (j=0;  j<mm2;  j++) 

{ 

if(mtfy[j]  !=  0.0)  if(py0  !=  0.0) 

py[j]=sqrt((re[j3*re[j]+imjj]*im[j])/py0)/mtfy[j]; 

} 

calc_ro(); 

initmtfQ; 

plot_mtf(); 

/*init_graph(); 

plot_graph();  *********  use  these  to  see  Isfs  */ 

} 

free(subframe); 

closeshutterO; 

} 

/* - Expose  CCD - */ 

void  expose(long  tm) 

{ 

clrccdQ; 
open_shutter(); 
long_delay(tm); 
closesh  utter  0; 

} 

/* - long  delay - */ 

void  long_delay(long  tm) 

{ 

delayparam  =  1000; 
while  (tm  >  1000) 

{ 

delay_ms(); 
tm  =  tm  -1000; 

> 

delay_param  =  (unsigned  short)tm; 
delay_ms(); 

} 

/* - calculate  the  line  spread  functions - */ 

void  calc_lsf(unsigned  int  far  *sub_frame) 

{ 
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unsigned  int  pixel; 
int  x,y; 

float  lsfxmax=0.0,lsfymax=0.0; 
for  (x=0;  x<165;  x++)  lsfx[x]=0.0; 
for  (y=0;  y<  1 65;  y++)  lsfy[y]=0.0; 
for  (x=size-1 +xstart;  x>=xstart;  x-) 

{ 

for  (y = size-1 +ystart;  y>=ystart;  y--) 

{ 

pixel =*(sub_frame++); 

Isfxfx]  =  lsfx[x]  +  (float)pixel; 

Isfyfy]  =  lsfy[y]  +  (float)pixel; 
if(lsfx[x]  >  Isfxmax)  lsfxmax=lsfx[x]; 
if(lsfy[y]  >  Isfymax)  lsfymax=lsfy[y]; 

} 

} 

lsfmax=lsfxmax; 

if(lsfxmax  <  Isfymax)  Isfmax  =  Isfymax; 

} 

/* - Subtract  off  the  background - */ 

void  background  () 

{ 

int  x.y.edge; 

float  sumx=0.0,sumy=0.0,ave_x1  ,ave_x2,ave_y1  ,ave_y2, aveback; 
edge=(int)(size/10.0);  /*this  uses  20%  to  find  background*/ 

/* — Find  the  average  background - */ 

for  (x=xstart;  x<xstart+edge;  x++)  sumx=sumx+lsfx[x]; 
ave_x1  =sumx/edge;  sumx=0.0; 

for  (x=xstart+ size-edge-1 ;  x<xstart+size;  x++)  sumx=sumx+lsfx[x]; 
ave_x2 = sumx/edge; 

for  (y=ystart;  y<ystart+edge;  y++)  sumy=sumy+lsfy[y]; 
ave_y1  =  sumy/edge;  sumy=0.0; 

for  (y=ystart+size-edge-1;  y<ystart+size;  y++)  sumy=sumy+lsfy[y] 
ave_y2 = sumy/edge ; 

aveback  =((ave_x1  +  ave_x2)  +  (ave_y1  +ave_y2))/4.0; 

/* — subtract  off  the  background - */ 

Isfmax  =  Isfmax  -  aveback; 

for  (x=xstart+size-1;  x>  =xstart;  x-) 

Isfxfx] = lsfx[x]-((x-xstart)  *(ave_x2-ave_x1  )/(size-edge) + ave_x1 ) ; 
for  (y=ystart+ size-1;  y>=ystart;  y--) 

Isfyfy] = Isfy[y]-((y-ystart)*(ave_y2-ave_y1  )/(size-edge) + ave_y1 ); 
/*  printf("\nThe  average  background  is  %f\n",aveback);*/ 
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> 

/********************pjp^  ppy******************************^ 

void  fft(int  m,  int  sign,  float  re[],  float  im[]) 

{ 

int  n,  j,  jl,  ndiv2,  nl,  n2,  k,  i,  ip,  npts; 
int  le,  leO,  lei,  I; 

float  pts,  t,  ure,  uim,  wre,  wim,  tre,  tim; 
double  ang,  pi; 

pi  =  4.0*atan(1.0); 

n  =  (int)(pow{2.0,  (double)m)  +  1 .0e-10); 

j  =  1; 

jl  =  0; 

nl  =  n  -  1; 

n2=  n  -  2; 

ndiv2  =  n/(int)2; 

for(i=0;  i<=n2;  i++) 

{ 

if(i  <  jl) 

{ 

t  =  re[j1];  /*  you  might  do  better  here  with  pointer  operations 

*/ 

re[j1]  =  re[i]; 
re[i]  =  t; 
t  =  im[j1J; 
im[jl]  =  im[i]; 
im[i]  =  t; 

} 

k  =  ndiv2; 
while(k  <  j) 

{ 

j  -=  k; 
k  /=  2; 

} 

j  +  =  k; 
ji=j-i; 

> 


for(l=1;  l<=m;  I++) 

{ 
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ure  =  (float)  1; 
uim  =  (float)O; 
ang  =  pi/(double)le; 
wre  =  (float)  cos  (ang); 
wim  =  (float)sin(ang); 
Ie0=le; 

Ie1=le-1; 
le  +=le; 


for(j=0;  j<=le1;  j  +  +) 

{ 

for(i=j;  i<=n1;  i+=le) 

{ 

ip  =  i  +  leO; 

tre  =  re[ip]*ure  -  im[ip]*uim; 
tim  =  re[ip]*uim  +  im[ip]*ure; 
re[ip]  =  re[i]  -  tre; 
im[ip]  =  im[i]  -  tim; 
re[i]  +=  tre; 
im[i]  +=  tim; 

} 

t  =  ure*wre  -  uim*wim; 
uim  =  ure*wim  +  uim*wre; 
ure  -  t; 


} 

} 

if(sign<0) 

{ 


pts  =  1/(float)n; 

for(i=0;  i<n;  i++) 

{ 

re[i]  *=  pts; 
im[i]  *=  pts; 

} 


> 

pQ**  *★*★★★*★**★*★  *★★***★★★★★*★*★★★  I 


void  calc_ro() 

{ 

int  j; 

float  z.zold.pj, slope; 


if  (px[1  ]  >  1.0)  if  (px[1]  <=  0.0)  printfHnX  DATA  POINT  BAD"); 
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else 

{ 

z  =  log  (-log(px[1  ])); 
for  (j=2;  j<mm2;  j  ++) 

{ 

pj  =  px[j]; 

zold  =  z; 
if  (pj<1.0) 

{ 

z  =  log  (-log  (pj)); 

slope  =  (z-zold)/log  ((float)  (j)/(float)(j-1)); 
if  (z  >  0.0)  goto  ro; 

} 

} 

} 

ro:  rox  =  rocal*exp(-zold/slope  +  log(j-1)); 

if  (py [1  ]  >  1.0)  if  (py [1  ]  <=  0.0)  printf("\nY  DATA  POINT  BAD"); 
else 

{ 

z  =  log  (-log(py [1  ])) ; 
for  (j  =  2;  j<mm2;  j  +  +  ) 

{ 

Pi  =  pyDl; 

zold  =  z; 
if  (pj<1.0) 

{ 

z  =  log(-log(pj)); 

slope  =  (z-zold)/log((float)  (j)/ (float)  (j-1 )) ; 
if  (z  >  0.0)  goto  ro2; 

} 

} 

} 

ro2:  roy  =  rocal*exp(-zold/slope  +  log  (j-1)); 

} 

/* . initialize  for  plotting--MTF's - */ 

void  init_mtf() 

{ 

int  TRUE=1 ; 
int  x,y; 

for  (x=0;  x<mm;  x+  +  ) 
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if(mtfx[x]  >  mtfmax)  mtfmax=mtfx[x]; 
for  (y=0;  y<mm;  y++) 

if(mtfy[y]  >  mtfmax)  mtfmax=mtfy[y]; 

setvideomode(  _VRES16C0L0R); 
setviewport(0,39,639,439); 

_setwindow(TRUE,0,0,mm+mm,  mtfmax); 
setcolor(15); 
moveto_w(0.,0.); 


/* - plot  modulation  transfer  functions - */ 

void  plot_mtf() 

{ 

int  x,y; 

for  (x=mm2;  x<mm;  x++) 

{ 

_lineto_w((double)(x-mm2),(double)mtfx[x]); 

} 

for  (x=0;  x<mm2;  x++) 

{ 

_lineto_w((double)(x+mm2),(double)mtfx[x]); 

} 

for  (y=mm2;  y<mm;  y++) 

{ 

_lineto_w((double)(mm2+  y),(double)mtfy[y]); 

} 

for  (y=0;  y<mm2;  y++) 

{ 

_lineto_w((double){mm+  mm2  +  y),(double)mtfy[yj); 

} 

_settextposition(1,1); 
printf("%s\n  %s”, TITLE, _ DATE _ ) ; 

_settextposition(4,1 ); 
printf("rox  =  %f",rox); 

_settextposition(4,38); 
printf("roy  =  %f",roy); 

_settextposition  (40, 30) ; 

_outtext(HHit  any  key  to  exit:"); 
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/*  getchO;  there  is  also  a  getch  in  subexpose,  only  use  one*/ 

} 


/* - initialize  for  plotting-lsf’s - */ 

void  init_graph() 

{ 

int  TRUE=1; 

_setvideomode{  _VRES16COLOR); 

_setviewport(0, 39, 639,439); 

_setwindow(TRUE,0,0,size+size,lsfmax); 

_setcolor(15); 

_moveto_w(0.,0.); 

} 

/* - plot  line  spread  functions - */ 

void  plot_graph() 

{ 

int  x,y; 

for  (x=xstart;  x<xstart+size;  x+  +  ) 

{ 

Jineto_w((double)(x-xstart),(double)lsfx[x]); 

} 

for  (y=ystart;  y<ystart+size;  y++) 

{ 

_lineto_w((double)(size+  y-ystart),(double)lsfy[y]); 

} 

_settextposition(1 ,1 ); 

printf ("%s\n  %s\n\n", TITLE, _DATE _ ) ; 

settextposition  (40, 30)  ; 

_outtext("Hit  any  key  to  exit:"); 

/*  getchO;  there  is  also  a  getch  in  subexpose,  only  use  one*/ 

} 
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APPENDIX  C  ARTIFICIAL  STAR  PROGRAM 


/*Uses  an  artificial  star  to  test  the  MTF.c  program*has  1st  plot  also**/ 
#define  TITLE  "Plot  an  artificial  star  MTF" 


#define  AUTHOR  "By  W.  J.  Rail" 

y********QujQj^  Q*************************************^ 


#include  <stdio.h> 

#include  <stdlib.h> 

#include  < graph. h> 

#include  <math.h> 

#include  <sys\types.h> 

int  xsize=165,  ysize=192,  xstart=0,  ystart=0; 
float  lsfx[165],lsfmax=0.0,lsfy[192]; 

int  sign=1,  m=8,  mm=2,  mm2; 

fl oat  mtfx [ 1 024] ,  mtf  max = 0.0, mtfy [ 1 024] ,  re [ 1 024] , im  [ 1 024] , px[51 2] , py [51 2] ; 
float  rocal  =  1 .0,rox,roy; 

void  make_star(); 
void  background  (); 

void  fft(int  m,  int  sign,  float  re[],  float  im[]); 
void  calc_ro(); 
void  init_mtf(); 
void  plot_mtf(); 

/*to  see  Isf’s  use  these  two  functions  instead  of  mtf  ones*/ 
void  init_graph(); 
void  plot_graph(); 


/* - main  routine - */ 

void  main() 

{ 

float  pxO.pyO; 
int  i,j,x,y; 

for  (i=0;  i<m-1;  i++)  mm  =  2*mm; 
mm2 = (int)  {(float)  (mm)/2.0) ; 
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make_star(); 
background  (); 


for  (x=xstart+xsize-1;  x>=xstart;  x--)  re[x]  =  lsfx[x]; 
sign,  re,  im); 

for  (x=0;  x<mm;  x++)  mtfx[x]  =  re[x]; 

/*find  the  x  power  spectrum***************/ 
pxO  =  re[0]*re[0]  +  im[0]*im[0]; 
for  (j=0;  j<mm2;  j++) 

{ 

if(mtfx[j]  !=  0.0  )  if(px0  !=  0.0) 

px[j] = sqrt((re  [j] *re  [j  J + im  [j]  *im[j])/pxO)/mtfx[j] ; 

} 

^reinitialize  the  FFT  arguements  before  finding  the  Y  FFT*/ 
for  (i=0;  i<1024;  i++) 

{ 

re[i]  =  0.0; 
im[i]  =  0.0; 

} 

for  (y=ystart+ysize-1 ;  y>=ystart;  y--)  re[y]  =  Isfy [y]; 
fft(m,  sign,  re,  im); 

for  (y=0;  y<mm;  y++)  mtfy[y]  =  re[y]; 

/*find  the  y  power  spectrum***************/ 
pyO  =  re[0]*re[0]  +  im[0]*im[0]; 
for  (j=0;  j<mm2;  j  +  +) 

{ 

if(mtfy[j]  !=  0.0)  if(py0  !=  0.0) 
py[j]=sqrt((re[j]*reO]+im[j]*im[j])/pyO)/mtfyIj]; 

} 

calc_ro(); 

init_mtf(); 

plot_mtf(); 

/*  init_graph(); 
plot_graph();  */ 

_setvideomode(  _TEXTC80); 

} 

/* - make  star - */ 

void  make_star() 
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int  x,y; 

float  pixel,lsfxmax=0.0,lsfymax=0.0; 
for  (x=xstart+xsize-1;  x>=xstart;  x--) 

{ 

for  (y=ystart+ysize-1 ;  y>=ystart;  y--) 

{ 

pixel=4095*exp(-0.1*(float)(((x-82)*(x-82)+(y-96)*(y-96)))); 

lsfx[x]  =  lsfx[x]  +  pixel; 

lsfy[y]  =  Isfy  [y]  +  pixel; 

if(lsfx[x]  >  Isfxmax)  lsfxmax=lsfx[x]; 

if(lsfy[y]  >  Isfymax)  lsfymax=lsfy[y]; 

> 

} 

lsfmax= Isfxmax; 

if(lsfxmax  <  Isfymax)  Isfmax  =  Isfymax; 

} 

/* - Subtract  off  the  background - */ 

void  background () 

{ 

int  x,y,edge; 

float  sumx=0.0,sumy=0.0,ave_x1  ,ave_x2,ave_y1  ,ave_y2,  aveback; 
edge=(int)(ysize/10.0);  /*this  uses  20%  to  find  background*/ 

/* — Find  the  average  background - */ 

for  (x=xstart;  x<xstart+edge;  x++)  sumx=sumx+lsfx[x]; 
ave  xl  =sumx/edge;  sumx=0.0; 

for(x=xstart+xsize-edge-1 ;  x<xstart+xsize;x++)sumx=sumx+lsfx[x]; 
ave_x2 = sumx/edge; 

for  (y=ystart;  y<ystart+edge;  y++)  sumy=sumy+lsfy[y]; 
ave_y1  =sumy/edge;  sumy=0.0; 

for  (y  =  ystart  +  ysize-edge-1 ;  y  <ystart  +  ysize;  y+  +  ) 
sumy = sumy + Isfy  [y] ; 
ave_y2= sumy/edge; 

aveback  =  ((ave_x1  +ave_x2)  +  (ave_y1  +ave_y2))/4.0; 

/* — subtract  off  the  background - */ 

Isfmax  =  Isfmax  -  aveback; 

for  (x=xstart+xsize-1;  x>=xstart;  x--) 

lsfx[x] =lsfx[x]-((x-xstart)*{ave_x2-ave_x1  )/(xsize-edge) + ave_x1 ) ; 
for  (y=ystart+ysize-1 ;  y>=ystart;  y--) 

Isfy[y]=isfy[y]-((y-ystart)*(ave_y2-ave_y1)/(ysize-edge)+ave_y1); 

/*  printf("\nThe  average  background  is  %f\n", aveback);*/ 
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/********************pjncj  ppy******************************/ 

void  fft(int  m,  int  sign,  float  re[],  float  im[]) 

{ 

int  n,  ],  jl,  ndiv2,  nl,  n2,  k,  i,  ip,  npts; 
int  le,  leO,  lei,  I; 

float  pts,  t,  ure,  uim,  wre,  wim,  tre,  tim; 
double  ang,  pi; 

pi  =  4.0*atan(1.0); 

n  =  (int)(pow(2.0,  (double)m)+1.0e-10); 

j  =  1; 

jl  =  0; 

nl  =  n  -  1; 

n2=  n  -  2; 

ndiv2  =  n/(int)2; 

for(i=0;  i<=n2;  i++) 

{ 

if(i  <  jl) 

{ 

t  =  re[j1  ] ;  /*  you  might  do  better  here  with  pointer  operations 

*/ 

re[j1]  =  re[i]; 
re[i]  =  t; 
t  =  im[j1  ]; 
im[j1]  =  im[i); 
im[i]  =  t; 

} 

k  =  ndiv2; 
while(k  <  j) 

{ 

j  -=  k; 
k  /=  2; 

} 

j  +=  k; 
ji=j-i; 

} 

le  =  1; 

for(l  =  1;  l<=m;  I++) 

{ 

ure  =  (float)  1; 
uim  =  (float)O; 
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ang  =  pi/(double)le; 
wre  =  (float)cos(ang); 
wim  =  (float)sin(ang); 
ieO=le; 
lei  =le-1 ; 
le  +=le; 


v 


( 


for(j=0;  j<=le1;  j++) 

{ 

for(i=j;  i<=n1;  i+=le) 

{ 

ip  =  i  +  leO; 

tre  =  re[ip]*ure  -  im[ip]*uim; 
tim  =  re[ip]*uim  +  im[ip]*ure; 
re[ip]  =  refi]  -  tre; 
im[ip]  =  im[i]  -  tim; 
re[i]  +=  tre; 
im[i]  +=  tim; 

} 

t  =  ure*wre  -  uim*wim; 
uim  =  ure*wim  +  uim*wre; 
ure  =  t; 


} 

} 

if(sign<0) 

{ 

pts  =  1/(float)n; 

for{i=0;  i<n;  i++) 

{ 

re[i]  *=  pts; 
im[i]  *=  pts; 

} 

} 

} 

!*** ************* ***Qg|0(j|g't0  pQ ***************************** *** i 

void  calc_ro() 

{ 

intj; 

float  z,zold,pj,slope; 


if  (ox[1]  >  1  0)  if  (px[1]  <=  0.0)  printfOnX  DATA  POINT  BAD"); 
else 
{ 
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z  =  log  {-log(px[1  ])); 
for  0=2;  j<mm2;  j  ++) 

{ 

pj  =  px[j]; 

zold  =  z; 
if  (pj<1.0) 

{ 

z  =  log(-log(pj)); 

slope  =  (z-zold)/log  ((float)  0)/(float)  0-1)); 
if  (z  >  0.0)  goto  ro; 

} 

} 

} 

ro:  rox  =  rocal*exp(-zold/slope  +  Iog0-1)); 

if  (py[1]  >  1.0)  if  (py [1  ]  <=  0.0)  printf('\nY  DATA  POINT  BAD"); 
else 

{ 

z  =  log  (-log(py[1  ])); 
for  0=2;  j<mm2;  j++) 

{ 

pj  =  py[j]; 

zold  =  z; 

if  (pj<1.0) 

{ 

z  =  log(-log(pj)); 

slope  =  (z-zold)/log  ((float)  0)/(float)  0-1)); 
if  (z  >  0.0)  goto  ro2; 

} 

} 

} 

ro2:  roy  =  rocal*exp(-zold/slope  +  Iog0-1)); 

} 

/* - initialize  for  plotting-MTF’s - */ 

void  init_mtf() 

{ 

int  TRUE=1; 
int  x,y; 

for  (x=0;  x<mm;  x++) 

if(mtfx[x]  >  mtfmax)  mtfmax=mtfx[x]; 
for  (y=0;  y<mm;  y++) 
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if(mtfy[y]  >  mtfmax)  mtfmax=mtfy[y]; 

setvideomode(  _VRES  1 6C0L0R) ; 
setviewport(0, 39,639,439); 

_setwindow(TRUE,  0,0,  mm +mm,  mtfmax); 
setcolor(15); 
moveto_w(0.,0.); 


/* - plot  modulation  transfer  functions - */ 

void  plot_mtf() 

{ 

int  x,y; 

for  (x=mm2;  x<mm;  x++) 

{ 

_lineto_w((double)(x-mm2),(double)mtfx[x]); 

} 

for  (x=0;  x<mm2;  x+  +  ) 

{ 

Jineto_w((double)(x+mm2),(double)mtfx[x]); 

} 

for  (y=mm2;  y<mm;  y++) 

{ 

_lineto_w((double)  (mm2  +  y),  (double)mtfy  [y]) ; 

} 

for  (y=0;  y<mm2;  y++) 

{ 

_lineto_w((double)(mm+  mm2  +  y),(double)mtfy[y]); 

} 

_settextposition(1,1); 

printf ("%s\n  %s", TITLE, _DATE_) ; 

_settextposition(4,1 ); 
printf  ("rox  =  %f",rox); 
settextposition  (4,38) ; 
printf  ("roy  =  %f",roy); 

_settextposition(40,30) ; 

_outtext("Hit  any  key  to  exit:"); 
getchO; 
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> 


/* - initialize  for  plotting — LSF’s - */ 

void  init_graph() 

{ 

int  TRUE=1; 

_setvideomode(  _VRES16COLOR); 

_setviewport(0, 39, 639,439) ; 

_setwindow(TRUE,0,0,xsize + ysize,  Isfmax) ; 

_setcolor(15); 

_moveto_w(0.,0.); 

} 

/* - plot  line  spread  functions - */ 

void  plot_graph() 

{ 

int  x,y; 

for  (x=xstart;  x<xstart+xsize;  x++) 

{ 

_lineto_w((double)(x-xstart),(double)lsfx[x]); 

} 

for  (y=ystart;  y<ystart+ ysize;  y++) 

{ 

Jineto_w((double)(xsize+  y-ystart),  (double)!sfy[y]); 

} 

_settextposition  (1 , 1 ) ; 
printf ("%s\n  %s\n\n", TITLE, _DATE_) ; 

_settextposition  (40, 30) ; 

_outtext("Hit  any  key  to  exit:"); 
getch(); 

} 
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APPENDIX  D  JITTER  PROGRAM 


/*jitter.c**********finds  the  cetroidjitter,  ro  and  f,  plots  isf******/ 

»  #define  TITLE  "Jitter  calculations  for  ro  and  f,  plots  Isf’s" 

#define  AUTHOR  "By  W.  J.  Rail" 

^*******QyjQ|^  Q**************************************^ 

/*#include  <alloc.h>  */ 

#include  <lynxx.h> 

#define  True  1 
#define  False  0 
#define  false  0 

#include  <stdio.h> 

#include  <stdlib.h> 

#include  <graph.h> 

#include  <math.h> 

#include  <time.h> 

#include  <sys\types.h> 

/*int  frame_ptr[FRAME_BYTES/2] ;  */ 
int  size=160,  xstart=0,  ystart=0; 
float  lsfx[165],lsfmax=0.0,lsfy[192]; 

float  xcsum =0.0,  xc2sum=0.0,  ycsum=0.0,  yc2sum=0.0;  /*for  jitter  calcs*/ 

int  n=1 ;  /*number  of  centroids  between  each  jitter  calc*/ 

float  dia=1.0,  focal=1000.0,  k=12.57E+6;/*wave  number,  500nm*/ 

float  xjit.yjit,  rocal  =  1.0,  rox,roy,fx,fy; 

void  init(); 

void  subexpose(long  tm,int  shutter); 
void  expose(long  tm); 
void  long_delay(long  tm)  ; 
void  calc Jsf (unsigned  int  far  *sub_frame); 

(  void  background{); 

void  centroid (); 
void  jitO; 
void  jit_rof{); 
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void  init_graph{); 
void  plot_graph(); 

/* - main  routine - */ 

void  main() 

{ 

initJynxxO; 

init(); 

cooler_off(); 

_setvideomode(  _TEXTC80); 

} 

/* - Initialize  program  info - */ 

void  init() 

{ 

int  cool  =  0,shutter=0; 
long  tm; 

printf ("Expose  time  in  ms.:"); 
scanf("%ld",&tm); 

printf("Cooler  0=OFF  1=ON:"); 
scant  ("%d",&cool); 
if  (cool)  cooler_on(); 

printf("Use  shutter  ?  (1  =Yes/0=No)"); 
scant ("%d",  &shutter); 

printf  ("Enter  size  of  subframe  in  pixels:"); 
scanf("%d",&size); 

printf("\nHit  any  key  to  quit:"); 

subexpose(tm, shutter) ; 

} 

/* - subExpose  CCD - */ 

void  subexpose  (long  tm.int  shutter) 

{ 

unsigned  int  *sub_frame; 
int  sub_frame_bytes,  i; 
double  deltime,time1,time2; 
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sub_frame_bytes  =  size*size  *  sizeof (unsigned  int); 
sub_frame  =  (unsigned  int*)malloc(sub_frame_bytes); 

if  (!  shutter)  open_shutter(); 

while  (  IkbhitO  ) 

{ 

for(i=0;  i<n;  i  +  +) 

{ 

if  (shutter)  expose(tm); 
else 
{ 

clrccd();  /*if  not  using  shutter*/ 
long_delay(tm); 

} 

x_param  =  (unsigned  char)  xstart; 
y_param  =  (unsigned  char)  ystart; 
size_param  =  (unsigned  char)  size; 
ptr_param  =  subframe; 
digitize_sub_frame(); 

calcjsf  (subframe) ; 
background  (); 
centroid  (); 

} 

jitO; 

jit_rof(); 

init_graph(); 

plot_graph(); 

> 

free(subframe); 
cl  osesh  utter  0; 

} 

/* - Expose  CCD - */ 

void  expose(long  tm) 

{ 

clrccd(); 

open_shutter(); 

long_delay(tm); 

close_shutterQ; 
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/* - long  delay - */ 

void  long_delay(long  tm) 

{ 

delay_param  =  1000; 
while  (tm  >  1000) 

{ 

delay  _ms(); 
tm  =  tm  -1000; 

} 

delay_param  =  (unsigned  short)tm; 
delay_ms(); 


/* - calculate  the  line  spread  functions - */ 

void  calc_lsf(unsigned  int  far  *sub_frame) 

{ 

unsigned  int  pixel; 
int  x,y; 

float  Isfxmax =0.0,  Isfymax =0.0; 
for  (x=0;  x<165;  x++)  lsfx[x]  =  0.0; 
for  (y=0;  y<165;  y++)  lsfy[y]=0.0; 
for  (x=size-1 -fxstart;  x>=xstart;  x--) 

{ 

for  (y = size-1 +ystart;  y>=ystart;  y-) 

{ 

pixel  =  *(sub_frame  +  +  ); 
lsfx[x]  =  lsfx[x]  +  (float)pixel; 
lsfy[y]  =  Isfy  [y]  +  (float)pixel; 
if(lsfx[x]  >  Isfxmax)  lsfxmax=lsfx[x]; 
if(lsfy[y]  >  Isfymax)  Isfy  max = Isfy  [y]; 

} 

} 

I  sf  m  ax = Isfxmax; 

if(lsfxmax  <  Isfymax)  Isfmax  =  Isfymax; 

} 

/* . Subtract  off  the  background — . */ 

void  backgroundQ 

{ 

int  x,y,edge; 

float  sumx=0.0,sumy  =  0.0,ave_x1,ave_x2,ave_y1,ave_y2,aveback; 
edge  =  (int) (size/1 0.0);  /’•'this  uses  20%  to  find  background*/ 

/* . Find  the  average  background . */ 
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for  (x=xstart;  x<xstart+edge;  x++)  sumx=sumx+lsfx[x]; 
ave_x1  =sumx/edge;  sumx=0.0; 

for  (x=xstart+size-edge-1;  x<xstart+size;  x++)  sumx=sumx+lsfx[x]; 
ave_x2 = sumx/edge; 

for  (y=ystart;  y<ystart+edge;  y++)  sumy=sumy+lsfy[y]; 
ave_y1=  sumy/edge;  sumy=0.0; 

for  (y=ystart+size-edge-1 ;  y<ystart+size;  y++)  sumy=sumy+lsfy[y]; 
ave_y2 = sumy/edge ; 

aveback  =  ((ave_x1  +  ave_x2) + (ave_y 1  +  ave_y2))/4.0; 

/* — subtract  off  the  background - */ 

Isfmax  =  Isfmax  -  aveback; 

for  (x=xstart+size-1;  x>=xstart;  x~) 

lsfx[x] = lsfx[x]-((x-xstart)  *(ave_x2-ave_x1  )/(size-edge) + ave_x1 ) ; 
for  (y=ystart+size-1;  y>=ystart;  y--) 

Isfy  [y] = Isfy  [y]  -  ( (y-ystart)  *  (ave_y2-ave_y  1  )/(size-edge) + ave_y  1 ) ; 

/*  printf("\nThe  average  background  is  %f\n",aveback);*/ 

} 

/* - Find  the  centroids - */ 

void  centroid () 

{ 

int  x,y; 

float  xcent,ycent,sumlx=0.0,sumly=0.0,sumxlx=0.0,sumyly=0.0; 


for  (x=xstart+ size-1;  x>=xstart;  x~) 

{ 

for  (y=ystart+size-1;  y>=ystart;  y--) 

{ 

sumxlx  =  sumxlx+  Isfx  [x]  *  (f  I  oat)  (x) ; 
sumlx  =  sumlx  +  lsfx[x]; 
sumyly  =  sumyly+  Isfy  [y  ]  *  (float)  (y) ; 
sumly  =  sumly  +  lsfy[y]; 

} 

} 

xcent=sumxlx/sumlx;  ycent=sumyly/sumly; 

/*printfHn  Centroid  is  xcent=%f  \n  ycent=%f,xcent,ycent);*/ 
/*  sum  centroids  for  jitter  calcs****************/ 
xcsum  =  xcsum  +  xcent;  xc2sum  =  xc2sum  +  xcent*xcent; 
ycsum  =  ycsum  +  ycent;  yc2sum  =  yc2sum  +  ycent*ycent; 


/****************** 


calculate  the  jitter 


A*************************/ 


void  jitQ 
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{ 

xjit  =  sqrt((xc2sum-xcsum)/(float)(n)); 
yjit  =  sqrt((yc2sum-ycsum)/(float)(n)); 
xcsum=0.0;  xc2sum=0.0;  ycsum=0.0;  yc2sum=0.0; 
printf("\n  Jitter  is  xjitter=%f  yjitter=%f, xjit, yjit); 

> 

/*******************ca|C(J|ate  ro  |  *from  t^e  jjttg,-****************/ 

void  jit_rof()  * 

{ 

float  a,b; 

a  =  xjit  *  k  /  focal; 

b  =  1.418  *  pow(dia, 0.333)  *  pow(a,2); 
rox  =  rocal  /  pow(b,0.6); 

a  =  yjit  *  k  /  focal; 

b  =  1.418  *  pow(dia, 0.333)  *  pow(a,2); 
roy  =  rocal  /  pow(b,0.6); 

fx  =  xjit; 
fy  =  yjit*roy; 

} 

/* - initialize  for  plotting . . */ 

void  init_graph() 

{ 

int  TRUE=1; 

_setvideomode(  _VRES16COLOR); 

_setviewport(0, 39, 639, 439) ; 

_setwindow(TRUE, 0,0, size + size.lsfmax) ; 

_setcolor(15); 

_moveto_w(0.,0.); 

} 

/* - plot  |jne  spread  functions - */ 

void  plot_graph() 

{ 

int  x,y;  ^ 

for  (x=xstart;  x<xstart+size;  x++) 

{ 

_lineto_w((double)(x-xstart),(double)lsfx[x]); 
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} 

for  (y=ystart;  y<ystart+size;  y++) 

{ 

_iineto_w((double)  (size  +  y-ystart) ,  (double)  Isfy  [y]) ; 

> 

_settextposition(1 ,1 ); 

'  printf("%s\n  %s\n\n", TITLE, DATE ) ; 

w-  _settextposition(4,1); 

printf("rox  =  %g\nfx  =  %g",rox,fx); 
settextposition  (4,38) ; 
printf("roy  =  %g",roy); 

_settextposition  (5,38) ; 
printfffy  =  %g",fy); 


settextposition  (40,30) ; 
outtextfHit  any  key  to  exit:"); 

/*  getchO;  subexpose  also  has  one,  only  use  one*/ 

} 


( 

4 
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APPENDIX  E  FRAME  RATE  MEASUREMENT  CODE 


/*time.c****times  digitize, calclsf  and  background,  cetroid, jitter*/ 

#define  TITLE  "Time  calculations" 

#define  AUTHOR  "By  W.  J.  Rail" 

/♦♦Quick  C  ******sample  of  time  code********************/ 
#include  <time.h> 

while  (  IkbhitO  ) 

{ 

timel  =  clockO;  /*  start  clock  */ 
for(i=0;  i<  128;  i++) 

{ 

digitize_sub_frame(); 
calc_lsf  (sub_frame) ; 
background  (); 
centroid  (); 

} 

jit(i); 


time2=clock();  /*  stop  clock  */ 
deltime  =  (time2-time1)/CLK_TCK; 

printf("\nTime  for  128  frames  =%f  sec",  deltime); 
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1.  Defense  Technical  Information  Center  2 
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Alexandria,  VA  22304-6145 
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Monterey,  CA  93943-5002 

3.  Prof.  Donald  L.  Walters  7 

Department  of  Physics  (Code  61  We) 

Naval  Postgraduate  School 
Monterey,  CA  93943-5004 

4.  Cpt.  Ann  Slavin  1 

Code  PLVWE 

Phillips  Laboratory 
Kirtland  Air  Force  Base 
Albuquerque,  NM  87117 
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